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ABSTRACT 

The galaxy intrinsic alignment causes the galaxy ellipticity-ellipticity power spectrum 
between two photometric redshifts to decrease faster with respect to the redshift sep- 
aration Az^ , for fixed mean redshift. This offers a valuable diagnosis on the intrinsic 
alignment. We show that the distinctive dependences of the GG, II and GI correlations 
on Az^ over the range \Az^\ < 0.2 can be understood robustly without strong as- 
sumptions on the intrinsic alignment. This allows us to measure the intrinsic alignment 
within each conventional photo-z bin of typical size > 0.2, through lensing tomography 
of photo-z bin size ~ 0.01. Both the statistical and systematical errors in the lensing 
cosmology can be reduced by this self-calibration technique. 
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1 INTRODUCTION 

The galaxy intrinsic alignment (lA) is one of the major sys- 
tematical errors of cosmic shear measurement. The intrinsic 
alignment of physically close galaxy pairs is correlated due 
to the tidal force arising from the correlated large scale 
structure and thus induces the intrinsic ellipticity-intrinsic 
ellipticity correlation (the II correlation). Due to the same 
reason, it is also correlated with the ambient matter distri- 
bution, which lenses background galaxies within sufficiently 
small angular separation and hence couples the shape of 
background galaxies with the shape of foreground galaxies, 
even if the galaxy pairs are widely separated in redshift. 
This induces the gravitation al shear-intrinsic ellipt icity cor- 
relation (the GI correlation. iHirata fc Seliakll200j 'l. Various 
methods have been propo sed to correct for the II correlation 
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dences can be understood robustly over the range jAz^l < 
0.2 without heavy lA modeling. This allows us to self- 
calibrate the intrinsic alignment with basically no assump- 
tions on the intrinsic alignment, through lensing tomogra- 
phy of fine photo-z bin size ~ 0.01 within each conventional 
photo-z bin of size > 0.2. With this self-calibration tech- 
nique, we no longer need to throw away the auto-correlation 
measurement of each conventional thick photo-z bins. It 
will bring in two-fold improvement on cosmology. First, it 
direc tly improves the c osmological constraints by 0(10%) 
(Ta kada fc Whitell2004 ) , since the auto-correlation measure- 
ment can now be safely included. The price to pay is orders 
of magnitude more correlation measurements due to many 
more photo-z bins. Second, it extracts valuable information 
on the intrinsic alignment within each thick photo-z bin, 
which helps calibrate the intrinsic alignment in the cross cor- 
relation meas urement. This new proposal, along with exist 
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Here we point out a new possibility to self-calibrate 
the intrinsic alignment. It is known that the lensing signal 
GG and the contaminations II and GI have different depen- 
dences on the pair separation Az^ . This has motivated the 
proposals to reduce the II correlation by using only cross 
correlation between th ick photo-z bins of size > 0.2 (e.g. 
iTakada fc Whitell2004l ). 

What we will show in this paper is that, these depen- 



ing ones (e .g. 'King fc SchneideiJ' 20o3: iBridl e fc Kind 
Joachimi fc Schneider 2008; Zhand 20081 : iBcrnsteiJ 
Joachimi fc Schnci derll2009l , I20 1Q) , demonstrate the rich in- 
formation brought in by the photo-z measurement to reduce 
systematical errors and to improve cosmological constraints, 
which is otherwise lost in the usual lensing tomography of 
thick photo-z bins. 

The current proposal is complementary to other model 
independe nt methods, especially the nulling technique pro- 
posed by IJoachimi fc Schneidej (2008, 2009, 2010). Our 
method takes advantage of the characteristic Az^ depen- 
dences of the lensing signal and the intrinsic alignment, for 
a fixed mean source redshift z'^ , to separate the two. Over 
the interesting range |Az^| < 0.2, since the lensing geom- 
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etry kernel is wide, the lensing signal barely changes as a 
function of Az^ while both II and GI change significantly. 
On the other hand, the nulling technique uses the charac- 
teristic dependence of the lensing signal, arising from the 
lensing geometry kernel, to suppress lA by a proper weight- 
ing function of z^ . It relies on significant variation of the 
lensing kernel with respect to over sufficiently large red- 
shift range. Hence the two methods are highly complemen- 
tary and we expect significant improvement in calibrating 
the intrinsic alignment by combining the two methods. 



2 THE AZ^ DEPENDENCES 

We work on the power spectrum C"''^ {£, Zi , Z2) between a 
property a at photo-z and another property /3 at Z2 ■ 
Here, the superscripts (and subscripts sometime) a, /3 — 
G, I, g. "G" denotes the lensing convergence k or the under- 
lying 3D matter distribution, "I" the E-mode intrinsic ellip- 
ticity and "g", the 2D or 3D galaxy number over-density. 
Throughout this paper, we focus on the Az^ = Z2 — zf 
dependence with fixed multipole £ and the mean redshift 
z^ = {zi + Z2)/2. Thus we will use the notation 

C''^{Az^\£,z- 



^aS I n P P\ 
C [l, Zi , 22 ) 



(1) 

and often neglect the arguments i and z^ . There is an 
obvious symmetry C^l^iAz^) = C'^°'{-Az^). It tells us 
dC^/dAz^lo = 0, a result which will become useful later. 
Unless otherwise specified, we wiU fix z^ = 1.0 and ^=10^ 
both are typical choices in weak lensing statistics. We adopt 
the Limber approximation to evaluate C"'^ in the standard 
flat ACDM cosmology. 
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Here, A\p{k,z) is the corresponding 3D power spectrum 
variance. x(.z) s-nd H[z) are the comoving angular diam- 
eter distance and the Hubble parameter respectively. The 
weighting function 
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) = W^{z,z^)Wp{z,z^) 



Wg{z,z'"} = H-'{z) I WL{z,z,)p{z,\z'')dz, 



(3) 
(4) 



Wi{z,z'') = W,(z,z^)=p{z\z'') . (5) 

Here, Wl{z, Zs) is the lensing kernel for a source at 2^ and a 
lens at z. p{z\z^) is th e photo-z PDF, model ed as the sum 
of two Gaussians fe..g. iMa fc Bernstein|[200^ . 
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Pent is the fraction of outlier galaxies, whose true redshift is 
biased by a factor /bias- We adopt 0-1,2(2^) = 0.05(1 -f z^) 
and /bias ~ 0.5. This toy model roughl y represents p(z|z^ ~ 
1) of a stage IV lensing survey (e.g. iBernstein k, Huterej 
|2009|'1. Stage IV da rk energy surveys require Pcat < 0.1% 
( Hearin et al]|2010l 'l. However, we adopt a much more con- 
servative Peat = 2%. Due to the possible photo-z scatters. 




Figure 1. The power spectra C°'l^{Az^\l = 10^,2^ = 1). We 
adopt the SB09 model to evaluate the II and GI correlations. Both 
II and GI vary with the pair separation Az^ in ways significantly 
different from the lensing signal, a key for the intrinsic alignment 
diagnosis. 



especially the catastrophic error, 

^1 J ^2) can be non- 
zero even if z[ < and C"^ {£, z[ , Z2) can be non-zero 
even if zC > Z2 ■ So the GI contamination that we re- 
fer throughout the paper is actually the sum of the two. 

To proceed, we adopt the lA model of 
ISchneider fc Bridlel (|2009l ') (hereafter SB09), based on 
the halo model prescription, as our fiducial model. SB09 
provides a fitting formula which allows us to conveniently 
perform the numerical calculation. We adopt the fiducial 
parameters of th e fitting formula listed in SB09 (and also 

r and hence C^^ and 



iKirk et al.ll2010l ') to calculate Ajj, 
C°'. 

C"^ show distinctive dependence on Az^ (Fig- HI - Both 
II and GI vary by ~ 20% from Az^ = 0.0 to Az^ = 0.1 
and by ~ 60% from Az^ = 0.0 to Az^ = 0.2. In sharp 
contrast, C'^'^ only decreases by 3% from Az^ = 0.0 to 
Az^ = 0.2. In another word, the Az^ dependences of the 
II and GI correlations are more than an order of magnitude 
stronger than that of the signal GG. This suggests that, if 
the intrinsic alignment contaminates the galaxy ellipticity- 
ellipticity power spectrum C''^' by a few percent, it would 
cause significantly different Az^ dependence, comparing to 
the case where the intrinsic alignment is ignored, in 

C^^\Az^) = C''^(A2^) (7) 
+ C"(Az^) + C'^'iAz'') + C"^{Az'') . 

The II term obviously causes C'^' to decease faster than the 
case without the II term. Since the GI term (= C^' + C"^) 
has a negative sign, the increment of its amplitude with 
Az^ means that it also causes C*-^' to decrease faster. Thus 
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Figure 2. Diagnose the intrinsic alignment through the A.z^ 
dependence of the galaxy ellipticity-ellipticity power spectrum. 
The data points result from the fiducial SB09 model and the 
error estimation is for LSST. The two dot lines have 50% weaker 
or stronger intrinsic alignment respectively. The solid line is the 
ideal case of no intrinsic alignment. The two dash lines are the 
results of an intrinsic alignment toy model with 7 = 1/2. The 
intrinsic alignment of the lower one is a factor of 2 stronger than 
the upper one. Results of other values of 7 have similar behavior. 
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Figure 3. The weighting function Wafiiz, /S.z^ = 1). The 
normalizations are fixed only within each o/3. The lines peak at 
left (black lines), middle (red lines) and right (blue) are Wqg^ 
Wgi + Wig S'li'i respectively. Solid, dot and dash lines have 

Az^ = 0.0,0.1,0.2 respectively. The small bumps at z ~ 0.5 of 
Wgi + Wig are caused by the combined effect of the photo-z 
outlier at 2 ~ 0.5 and the lensing peak at 2: ~ 0.5. They have 
only weak dependence on Az^ . 



the intrinsic alignment always causes C^^' to decrease faster 
with respect to Az^ (Fig. [2)|. For example, for the SB09 
model, C^^' decreases by 10% to Az^ = 0.2, comparing to 
the 3% decrease in C^^^ . 

To assess the generality of the above behavior, we also 
investigate a toy model. In this toy model, the intrinsic 
alignment has a bias (with respect to matter distribution) 
bim{k,z) (X [l + Al,{k,z)y ( 7 G [0,1/2]) and the cross 
correlation coefficient r — —1. Here, A'^{k.z) = AQQ{k,z) 
is the 3D matter power spectrum. For this set up, Ajj = 
6|„A^ and Afg = Aqj = fc/^rA^ = -bim^m- There is 
no solid physics behind this toy model. But roughly speak- 
ing, 7 = mimics a class of lA models with linear depen- 
dence on the matter overdensity and 7 = 1/2 with quadratic 
dependence. By choosing different 7 we can cover a wide 
range of lA scale and redshift dependence. Over 7 G [0, 1/2], 
we confirm the behavior that the intrinsic alignment causes 
C^^' to decrease faster with respect to Az'^ (Fig. [SI- 

This faster than usual decrease is a smoking gun of the 
intrinsic alignment. However, this smoking gun survives only 
when the photo-z performance is reasonably good. For exam- 
ple, if the photo-z measurement is completely wrong, with 
no correlation with the true redshift at all, the Az^ depen- 
dence would vanish. Fortunately, Fig.[2]shows that, for typi- 
cal photo-z performance accessible to stage IV projects, the 
Az^ dependence is largely preserved, even at Az^ < 0.1. 
The usual lensing tomography with coarse redshift bins of 



size > 0.2 thus misses the valuable information of the Az^ 
dependence at Az^ < 0.2. Such information can be recov- 
ered by lensing tomography with fine bin size ~ 0.01 within 
each coarse bin, albeit requiring two orders of magnitude 
more computations. 

LSST and other surveys of comparable capabil- 
ity are able to measure the Az^ dependence over 
Az^l < 0.2 to desired accuracy. To a good approxi- 
mation, C^^ at different Az^ samples the same cosmic 
volume (Fig. |3J and thus shares the same cosmic vari- 
ances. The relative differences are then ^ {£A£fsky)~^^^ ~ 
0.3%(10V^)^''^(100/A.^)i/VsVy . not a limiting factor to 
measure the Az^ dependence. 

The random galaxy shape shot noise dominates over the 
cosmic variance at ^ > lO'^, for a photo-z bin of size 0.01 at 
z^ = 1 wit h 2.5 X 10^ galaxies, a typical number for LSST 
(jZhan et a l. 2009). Furthermore, since shot noises are uncor- 
related at different Az^ , it becomes more difficult to mea- 
sure the Az^ dependence. However, since the Az^ depen- 
dence has only weak dependences on £ and , we can adopt 
a wide bin size A£ = 500 and average over 10 of these mea- 
surements of the same Az'^ but diflferent z'' G [0.95, 1.05] to 
beat down the shot noise to percent level, sufficient to diag- 
nose lA and distinguish between some lA models of interest 
(Fig.©. 
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3 UNDERSTANDING THE AZ^ 
DEPENDENCES 

Eq.[2]shows that, C^^^l, Az^ , z^) of the same i and sam- 
ples the same 3D clustering A^^, but with different weight- 
ing function Wap, which is the only function in the integrand 
depending on Az^ . Hence the Az^ dependence in C"^ , es- 
pecially the ratio C{Az )/C{Az = 0), should be mainly 
determined by the Az^ dependence in Wap, but not A^^. 
This simple fact turns out to be highly valuable. (1) It al- 
lows us to understand the Az^ dependence in C'^'^ to 0.1% 
accuracy f i]3.1|) . (2) It allows us to understand the Az^ de- 
pendences in C'' (3321) S'lid C'^' ( t|3.3[) without strong as- 
sumptions on the intrinsic alignment, since Wap does not 
depend on the property of the intrinsic alignment. These 
are keys towards a model-independent intrinsic alignment 
self-calibration (Q. 



3.1 The Ajz^ dependence in C"^ 



Wgg{z,Az' 



peaks at half the distance to the source. 



Although the peak amplitude is sensitive to z , it only 
weakly depends on Az^ and only decreases by a few per- 
cent from Az^ = to Az^ = 0.2. This explains the weak 
dependence of C'^'^ on Az^ , which can be well described by 
the Taylor expansion around Az^ = up to second order, 
since 9C°''/aA2^|o = 0, 



C°c:[AzP =0) 



fGG{l,z''){Az''f 



(8) 



a"C^"(A^^)/a(Az^)^|o 
CGG(Az-P = 0) 



We find that, /gg £ (0.5, 1.0) for t G [20, 4000] and z^ = 1. 
At ^ = 10^, /gg = 0.7. Eq. [8]is accuracy to 0.1% up to 



Az^ 



0.3. 



3.2 The Az^ dependence in the II correlation 

Wii{z, Az^ , z^) peaks sharply at a true redshift Zpcak (in 
our case, Zpeak — 1). The peak position is insensitive to Az^ , 
although it is sensitive to z^ , as long as the photo-z measure- 
ment is sufficiently accurate. The peak amplitude decreases 
with jAz^l quickly, causing the sharp decrease in C^' at 
jAz^l ^ 0.1. This behavior is well known in the literature 
and has been applied to reduce the II correlation by using 
only the cross power spectra between thick photo-z bins, 
corresponding to the limit of |A2;^| ^0.1. Instead, here we 
explore the II behavior at |A2^| < 0.2, namely within each 
thick photo-z bin of conventional lensing tomography. 

Since Wii has a dominant and narrow peak whose po- 
sition Zpcak is insensitive to Az^ , 




Figure 4. The accuracy of the scaUng relation EQ.l8l [Tolllll fc [T2l 
The two lines labeled with "GG" are the left and right hand sides 
of Eq. [8] Eq.|8]is accurate to 0.1% and thus the two lines overlap. 
The intrinsic alignment models demonstrated are the SB09 model 
and the toy model with 7 = 1/2. The toy models with 7 < 1/2 
result in better accuracy of the scaling relation. Eq. 1101 &: llll are 
accurate to ~ 1% so the lines labeled with "gg,ll,lg" overlap. The 
upper solid line (blue) is C'^^ + C^'^ . The two dash (magenta) 
lines labeled with "GI+IG" are the SB09 model (lower one) and 
the toy model with 7=1/2 (upper one). Eq. 1121 is accurate to 
better than ~ 5% at Az^ ^ 0.2. 



Furthermore, since Wn 
relation 



Wgg, we suggest the following 



C"{Az' 



) ~ Aii{£,z^)C<''{Az'^\£,z^) 



(10) 



meaning the same Az^ dependence in C" and C"". This 
relation is not only more accurate in general, but also more 
useful. The same lensing survey measures C^^{Az^) and 



thus tells the Az^ dependence in C^^, 



without external 



knowledge on p{z\z ) nor the intrinsic alignment. For ba- 
sically the same reason, we also have 



c'''{Az''\e, z'') ~ Aig{i, z'')C{Az''\e, 2-^) 



(11) 



Eq. [To] & \TT\ are accurate to 1%, for the SB09 model and 
the toy model (Fig. Q. To avoid modeling uncertainty in 
All and Ajg (and Agi in next subsection), we treat them 
as free parameters in the proposed self-calibration. 



C"(A2^|£,2^) 



27r^X(2peak)iy (apeak) .2 



Aij{k,z) 



Wii{z, Az^ , z^)dz 



(9) 



Here, Ajj(fc,z) is evaluated at fc = i/x{z) and z = ^pcak. 
The Az^ dependence is completely described by the last 
integral. In the limit that pcat ^ and the dependence of cti 
on z^ is negligible, the predicted dependence becomes exact. 



3.3 The Az^ dependence in the GI correlation 

The corresponding weighting function for the GI correlation 
iC^' + C"^) is Wgi + Wig- It peaks at redshift lower than 
z^ . When Az^ increases, the true redshift separation be- 
tween source and lens increases and Wgi + Wig increases. 
Thus the amplitude of the GI correlation {\C'^' + C"^\) in- 
creases with Az'^ . Due to the peak feature in Wgi + Wig, 
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we postulate 

C'=^^(Az^i£, z'') + C^^(Az-^|^, z"") ~ Aaiit, z^) 

X [c''»(A2-^|£,z-^) + C»'^(A2-^|^,2-^)] . (12) 

Since the peak is not as sharp as the one in Wn and the 
peak position does move with respect to Az^, and since the 
contribution from those outlier galaxies at 2 ~ 0.5 becomes 
non- negligible (Fig. [3]), the accuracy of Eq.[T2]is not as good 
as Eq. IIOI however, it still reaches ~ 10% up to t\z^ = 0.2 
(Fig. |4]). We thus take a conservative approach by restricting 
the self-calibration technique to Az^ < 0.2. For this reason, 
performing fine binning of size 0.01 within coarse bin of size 
> 0.2 suffices. However, if we can improve the accuracy of 
Eq. HH beyond Az^ = 0.2, binning the whole redshift range 
into fine bins of Az^ ~ 0.01 would be beneficial. 

The above results demonstrate that the scaling rela- 
tions [TOl [11] & [12] are indeed the manifestations of their cor- 
responding weighting functions, which do not rely on the 
lA properties. The scale and redshift dependences of the lA 
clustering property are significantly different between the 
toy model and the SB09 model, further suggesting the gen- 
erality of the above Aa^ dependences. Improvements over 
llOUlll fc ll2l if necessary, can in principle be achiev ed by the 
theory of Gaussian quadratures (|Press et al.lll997l ). 



with the presence of other errors, such as the PSF, which 
should have different t\z^ dependences. 

The intrinsic alignment may also induce a non- 
negligible B-mode shape distortion (e.g. iHevmans et al.l 
200i|) , whose power spectrum should have virtually the same 
t\z^ dependence as C®® and differ significantly from those of 
other B-mode sources. This offers an independent way to di- 
agnose lA from the B-mode shape distortion measurement, 
although extra modeling is required to use this B-mode mea- 
surement to correct for the E-mode I A. 

The self-calibration proposal rely on no external mea- 
surement nor strong lA assumptions. In this letter we 
present a concept study on its feasibility. Quantitative anal- 
ysis on its performance, along with comprehensive investi- 
gation on various complexities in realistic surveys, shall be 
carried out to robustly evaluate this proposal. 
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4 SELF-CALIBRATING THE INTRINSIC 
ALIGNMENT 

The discovered scaling relations (Eq. O \W\ \TT\ & [T^ can 
be conveniently plugged into existing framework of lens- 
ing tomography analysis to self-calibrate the intrinsic align- 
ment. Basically, lensing surveys allow for the measure- 
ment of the Az^ dependences in C*^' (Eq. [7]), C^^' (the 
ellipticity-galaxy density power spectrum) and C'^' (the 
galaxy density-galaxy density correlation pow er spectrum) 
(hu fc Jainll20^ : lBernstein|[2009l : IZhanell2008l ') Q 

C'-^'>{Az'') = C'^»(A2-^) + C"''^(A2^) + 2C^»(A2^) , 
C<^'(A2'^) = C'^iAz'') . (13) 

One can show that, with the aid of Eg. [8l [TOl [TTI fc [T2l 

measurements at 4 or more Az^ allow for simultaneous re- 
construction of GG, II, GI and free parameter fcG, Aqi, 
All and Aig. 

The total contamination (II+GI) in C*^' (Eq. [7|) can 
be measured with higher accuracy than II or GI, since both 
causes C'^' to decrease faster and are thus partly degener- 
ate. Basically, a 3% lA contamination in C'^' would double 
the Az^ dependence, observable by LSST (Fig. [2]) or other 
surveys with comparable capability. 

Although we focus on ^ = 10'^ and z^ = 1, the self- 
calibration proposal is applicable to other redshifts and an- 
gular scales. For example, it may work better at larger scales 
I ~ 10^, where shot noise is less an issue. It may also work 



^ The magnification bias can add non-negligible corrections pro- 
portional to 0*^9, C'^ & C^"^ to Eq.[T3l Since none of them is 
new unknown, it does not invalidate the self-calibration, although 
measurement error in the galaxy luminosity function does bring 
new uncertainties. 
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